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ABSTRACT 

In a previous paper (Ouyed et al. 2004) we presented a new model for soft gamma-ray 
repeaters (SGR), based on the onset of colour superconductivity in quark stars. In this 
. model, the bursts result from the reorganization of the exterior magnetic field following 

\^ , the formation of vortices that confine the internal magnetic field (the Meissner effect). 

\ Here we extend the model by presenting full 3-dimensional simulations of the evolution 

ly-^ \ of the inclined exterior magnetic field immediately following vortex formation. The 

' simulations capture the violent reconnection events in the entangled surface magnetic 

i-S^ ■ field as it evolves into a smooth, more stable, configuration which consists of a dipole field 

I ' aligned with the star's rotation axis. The total magnetic energy dissipated in this process 

o 



is found to be of the order of 10 erg and, if it is emitted as synchrotron radiation, peaks 



c/3 . typically at 280 keV. The intensity decays temporally in a way resembling SGRs and 

^ ' AXPs (anomalous X-ray pulsars), with a tail lasting from a few to a few hundred times 

the rotation period of the star, depending on the initial inclination between the rotation 
^ ' and dipole axis. One of the obvious consequences of our model's final state (aligned 

' rotator) is the suppression of radio-emission in SGRs and AXPs following their bursting 

era. We suggest that magnetar-like magnetic field strength alone cannot be responsible 
for the properties of SGRs and AXPs, while a quark star entering the "Meissner phase" 
is compatible with the observational facts. We compare our model to observations and 
highlight our predictions. 

Subject headings: gamma rays: bursts — X-rays: stars — stars: magnetic fields — 
stars: neutron — stars: quark star 



1. Introduction 



Soft 7-ray repeaters (SGRs) are sources of recurrent, short {t ~ 0.1 s), intense (L ~ 10^-lO^LEdd) 
bursts of 7-ray emission with a soft energy spectrum. The normal pattern of SGR activity are in- 
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tense activity periods which can last weeks or months, separated by quiescent phases lasting years 
or decades. The five known SGRs are located in our Galaxy or, in the case of SGR 0526-66, in 
the Large Magellanic Cloud. The two most intense SGR bursts ever recorded were the 5 March 
1979 giant flare of SGR 0526-66 (Mazets et al. 1979) and the similar 27 August 1998 giant flare 
of SGR 1900-1-14. The peak luminosities of these events (~ lO^-lO^LE^d) exceeded the peak lu- 
minosities of "normal" SGR bursts by a factor > 10^. Several SGRs have been found to be X-ray 
pulsars with an unusually high spin-down rate oi P/P ~ lO"^'^ s~^, usually attributed to magnetic 
braking caused by a super-strong magnetic field B > lO^^G, which leads to the interpretation that 
SGRs are magnetars (Golenctskij ct al. 1979; Duncan & Thompson 1992, Kouveliotou et al. 1998, 
Kouveliotou et al. 1999). In the magnetar model, the magnetic field is the likely provider of the 
burst energy. A common scenario assumes that magnetic stresses create a quake in the crust of the 
neutron star, which then ejects hot plasma Alfven waves through its rigid magnetosphere (Thomp- 
son & Duncan 1995; 1996). The magnetic field of such a star would have grown to magnetar-scale 
strengths because of strong convection during the collapse of the proto-neutron star core (Duncan 
& Thompson 1992; Thompson & Duncan 1993). 

1.1. Open issues in the magnetar model of SGRs 

In the magnetar model of SGRs, which is also that of Anomalous X-ray Pulsars (AXPs), the 
X-rays are ultimately powered by an internally decaying very strong magnetic field. However there 
are still a few open questions which in our opinion leave room for new models to be explored: 

• Despite numerous attempts, no magnetars have been detected at radio frequencies^ (Kriss et 
al. 1985; Coe et al. 1994; Lorimcr et al. 1998; Gaensler et al. 2001). It has been suggested 
that QED processes at high B, such as photon splitting, may preclude the electron/positron 
cascades necessary to produce radio emission (Baring&Harding 2001), or that pair production 
ceases above some critical magnetic field (Zhang & Harding 2000). These ideas or alternatives 
remain to be confirmed and are still debatable. 

• One might expect high-S radio pulsars to be more X-ray bright than low-i? sources, and 
to possibly exhibit AXP-like burst emission. However X-ray observations of 5 high-i? radio 
pulsars reveal luminosities much smaller (by a few orders of magnitudes) than those of AXPs 
(Pivovaroff et al. 2000; Gonzalez &: Safi-Harb 2003; McLaughlin et al. 2003; Gonzalez et 
al. 2004; Kaspi & McLaughlin 2005). This has lead to suggestions that high-S radio pulsars 
may one day emit transient AXP-like emission, and conversely that the transient AXPs might 
eventually exhibit radio pulsations (Kaspi &; McLaughlin 2005) - a notion yet to be confirmed. 

• Hints of massive (> 30-40 Mq) progenitors associated with AXPs and SGRs by recent ob- 



^Detecting radio pulsations may be difficult, given the small polar caps implied by the long spin periods. 
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servations (Gaensler et al. 2005) has led to the suggestion that pulsars and SGRs differ in 

their progenitor masses. It is also suggested that massive progenitors could lead to neutron 
stars with millisecond periods (Hcger ct al. 2004) which would comply with the magnetar 
model for SGRs (Duncan & Thompson 1992). This, however, leaves open the question of why 
high-B pulsars, formed from less massive progenitors (presumably with periods > 15 ms), 
possess magnetar-like field strengths. 

• All SGRs and AXPs known to date have spin periods between 6 and 12 seconds (Kaspi 2004). 
This clustering in period remains to be explained by the magnetar model. 

Here we explore an alternative model first presented by Ouyed et al. (2004) where one assumes 
that AXPs and SGRs are quark stars, rather than magnetars. While quark stars have yet to 
be found in nature, formation scenarios have been suggested in the literature (see Sect. 8.4 in 
Ouyed et al. 2004 and references therein). The qualitative idea is that the core of a neutron 
star reaches deconfinement densities, eventually leading to the conversion of the entire star to a 
quark star. In principle, the transition from hadronic matter into quark matter in the core of 
a neutron star can happen immediately during or after the supernova explosion, but also much 
later than that. Such a transition could occur in a smooth stable manner (e.g. Bombaci & Datta 
2000 and references therein) or in an explosive manner termed the "Quark-Nova" (Ouyed et al. 
2002; Keranen, Ouyed, &: Jaikumar 2005). In these formation scenarios, if the quark star is born 
shortly after the supernova, the emitted energy /radiation would be absorbed by the still expanding 
supernova ejecta, and there would be no detectable signature of the transition. The hadron-quark 
transition may also happen much later than the supernova explosion and could be induced, for 
example, by accretion if the neutron star is a member of a binary, or via spin-down in the case of 
isolated neutron stars. The new idea of this model, first presented in Ouyed et al. (2004), is that 
the quark star enters a superconductive phase, and subsequently experiences a "Mcissner phase" 
that triggers the reorganization of the star's magnetic field. Before going into more details we first 
describe properties of quark matter and the concept of color superconductivity. 



1.2. Superconductivity and Meissner effect in quark stars 

The discovery of asymptotic freedom, leading to the formulation of quantum chromodynamics 
(QCD) as the theory of strong interactions, was soon followed by the suggestion that matter at 
sufficiently high densities consists of a deconfined phase of quarks (Collins & Perry 1975). Only 
shortly afterwards it was pointed out (Barrois 1977) that the true ground state of cold dense 
quark matter exhibits color superconductivity (CSC), characterized by diquark condensation with 
an estimated energy gap A ~ 1 MeV between the highest occupied and the unperturbed quark 
state at the Fermi surface. Since this magnitude of the gap is rather small for phenomenological 
applications, CSC subsequently received little attention. The situation changed when reinvestiga- 
tions (Alford, Rajagopal, & Wilczek 1998; Rapp et al. 1998), using nonperturbative forces {e.g., 
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instanton- induced) , showed that the gap can be substantially larger, A ~ 100 MeV for moderate 

quark chemical potentials, fig ~ 350 MeV. Similarly large values are obtained from estimates based 
on perturbative calculations at asymptotically high densities (Pisarski & Risclikc 1999; Son 1999). 
Thus, from the practical point of view, the existence of color superconductivity in compact stars 
has (re-)emerged as an exciting possibility. 

The detailed properties of CSC matter relevant to astrophysical applications depend on the 
interplay of the quark chemical potential, the q-q interaction strength, and the bare masses of the 
(light) quarks u, d and s. In particular, for /Xg below the (constituent) strange quark mass, only 
u and d quarks are subject to Cooper pairing. The corresponding phase is known as 2-flavor CSC 
(2SC). In the idealized case where the quark chemical potential is much larger than the strange 
quark mass (rris), the latter becomes negligible and all three flavors exhibit likewise pairing. The 
preferred symmetry (breaking) pattern in this phase corresponds to the so-called color-flavor locking 
(CFL; Alford, Rajagopal, & Wilczek 1999), since the underlying diquark condensate is invariant 
only under simultaneous color and flavor transformations. In the present work, we will focus on 
the CFL phase (for a recent review and a more exhaustive list of references, cf. Schafer 2003). 
Associated with CSC is the critical temperature, fceTc ~ 0.5 A, above which pairing is washed out. 

One of the most interesting properties of an ordinary superconductor is the Meissner effect, 
i.e., the expulsion of magnetic flux from the superconductor (Meissner &: Ochsenfeld 1933). In the 
CFL phase, the gauge bosons connected with the broken generators obtain masses, which indicates 
the Meissner screening effect (Rischke 2000a,b). This is the heart of our model which we describe 
next. 

1.3. Our model 

Assume a quark star is born with a temperature T > and enters the CFL phase as it cools by 
neutrino emission (Keranen, Ouyed &; Jaikumar 2005). The CFL front quickly expands to the entire 
star followed by the formation of rotationally induced vortices, analogous to rotating superfluid ^He 

(the vortex lines are parallel to the rotation axis; Tilley&Tilley 1990). Via the Meissner effect, the 
magnetic field is partially screened from the regions outside the vortex cores. The system now 
consists of alternating regions of superconducting material with a screened magnetic field, and the 
vortices where most of the magnetic field resides. As discussed in Ouyed et al. (2004), this has 
interesting consequences on how the surface magnetic field can adjust to the interior field which is 
pinned to the vortices. Figure 1 shows the starting point for this reorganization of the magnetic 
field, just after the Meissner effect has fully aligned the magnetic flux inside the star with the 
rotation axis. Within a transition region, the magnetic field switches from the vertical interior field 
to that of an inclined dipole outside. A conservative assumption (from the energetics point of view) 
is that the field in the transition region is a potential field (minimum-energy configuration). 

In order to capture the complex non-linear dynamics of the system we need the help of 3-D 
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simulations. These simulations will allow us to follow the evolution and re-organization of the 
magnetic field induced by the vortices. 
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Fig. 1. — "Mock magnetic field lines" at t = in the xz-plane. The rotation axis is in the z direction 
inclined by an angle 9 with respect to the dipole axis. Initial state after the stars interior converts 
to the CFL phase thus causing the magnetic field to be confined inside the vortices. The inner 
circle represents the surface of the star, and between that and the outer circle is the region where 
surface currents cause the outer dipole field to adjust to the vertex-aligned inner field. Shown are 
isolines of the vector potential component Ay, which for our setup trace most features of magnetic 
field lines. 

The paper is presented as follows: In § 2 we describe and analyze the basic setup of our 
simulations. We calculate the synchrotron light curves that result from the simulations wherein we 
find remarkable similarities to SGR light curves. Mechanisms for the observed subsequent bursts 
and the quiescent phase in SGRs are then briefly discussed in terms of our model in § 3. In § 4 we 
discuss the model predictions and how it can account, at least in its current stage, for the points 
listed in § 1.1. Here we suggest a list of observations that could test our model. We conclude in 
S 5. 
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2. Simulations 

Two fundamental parameters of our model are the inclination angle 9 of the external dipole field 
relative to the rotation axis, and the initial ratio of gas to magnetic pressure at the surface. The 
simulations start with the interior magnetic field confined in vertical vortices for r < Rqs , but with 

a still unperturbed inclined dipole field for r > 1.3 /Zqs- The transition region Rqs < r < 1.3i?Qs is 
filled by a potential field and is bounded by two current layers (sec Fig. 1). The whole configuration 
was chosen by minimizing the total magnetic energy (see Appendix B). 

We solve the following set of non-ideal magnetohydro dynamical (MHD) equations, using the 
Pencil-Code (see e.g. Dobler et al. 2006)^. 

Ding _ _^ 

= -V-u (1) 



Dt 

Du 9^ f s , \ „^ JxB 



+ u ^'^u + ^VV ■ U + 2S -Vlng 

dA 

— = uxB- rinoj (3) 
Ds 

qT— = V ■{KVT)+'niiof +2gvS'' . (4) 

Here g, u, A, and s are density, velocity, magnetic vector potential, and specific entropy, respectively. 
Parameters and functions kept constant were Cp, <5grav5 and rj corresponding to specific heat, 
gravity potential, kinematic viscosity and magnetic diffusivity respectively. The remaining variables 
Cs, j, B, S and T represent the sound speed, electric current density, magnetic flux density, traceless 
rate-of-strain tensor, and temperature, respectively. 

The length scale is in units of the radius of the quark star, Rqs, density in units of the star's 
surface density, po, and time is in units of the spin-period, This implies that velocities are 

in units of Rqs^, while the magnetic field is in units of ./pqRqs^- The strength of the (dipole) 
magnetic field at the surface can be estimated for a given /3o as, B^ = SttPo/Pq, where Pq is the 
pressure at the surface of the star. Using hydrostatic balance and the perfect gas law this becomes, 

^0 ~ 7^= — X (5) 

vPo 

PO \ / 10km\ /Mqs X 1/2 

,106g/cm3j V ^Qs ; \Mq) 

In the equation above, po is the average density of the gas close to the surface of the star. This 
corona is supplied by fall-back material following the formation of the QS (Keranen, Ouyed, &: 



^http://www.nordita.dk/software/pencil-code - The Pencil Code is a high-order finite-difference code for solving 
the compressible hydromagnetic equations. 
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Jaikumar 2005) similar to what has been suggested in the supernova case (ChevaHer 1989). The 
density of the fallback matter, representative of the crust material of the parent neutron star, is 
estimated to be of the order of 10^ g cm~^. 

Note that our MHD equations (1)— (4) are non-relativistic. While near the surface of the quark 
star some aspects of the physics will be considerably changed by rclativistic effects, we expect that 
the overall dynamics will not be vastly different in a more realistic calculation. 

All figures shown here are for resolution 128^. Wc have run simulations at higher resolution, 
but found that they differ little as far as energetics and evolution are concerned. 

2.1. Evolution and reorganization of the surface magnetic field 

Figures 2-3 show the evolution of the exterior magnetic field as it adjusts to the vortex-confined 



The complicated structure of the surface magnetic field is clearly seen in the animations driven 
by the frequent magnetic reconnections as the surface field tries to align itself with the interior one 
(rotation axis). These random reconnection events would bear many similarities to the initial events 
{i.e. those near t = 0), but we expect them to be less energetic as the magnetic field slowly decays 
and weakens. Eventually, the magnetic field evolves into a stable configuration (see Fig. 2) after 
which the star enters a quiescent phase. 

The restructuring of the field in the transition region leads to an approximately spherical 
Alfven wave traveling outwards (see Fig. 2 for t = 1) and is more prominent in simulations with 
a stronger magnetic field (i.e. smaller (3). As the wave travels outwards it amplifies the magnetic 
field in certain regions causing them to undergo reconnection, which both distorts the wave and 
eventually damps it out. Furthermore, the regions that underwent reconnection appear to show 
slow oscillatory motions between the reconnection site and the surrounding gas ("breathing"). This 
can be seen in the series of diminishing pulses in Fig. 6 and the frequency of the pulses remains 
nearly constant (see Fig. 7). We note that these pulses appear more prominently in simulations 
with lower /3 and do not arise in simulations with P > 1. 



The magnetic energy released in the organization is shown in Fig. 6 and can be cast into a 
simple equation. 



field^. 



2.2. Energetics and emission 





^See the" Animations" link at http://www.capca.ucalgary.ca for movies of the simulations. 
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where a is the fraction of the surface magnetic field that decayed via reconnection events. Prom 
Fig. 4 we can infer the decay in magnetic energy to be roughly a = 0.4 - 0.6. 

Furthermore, if we assume the electrons to be co-rotating (7e ~ 1), then the intensity emitted 
by synchrotron processes from the simulated region goes as Is oc UeB"^. Also, assuming peak 
emission is at Up = 0.29z^c5 where Uc is the critical frequency (see Sects. 6.2 and 6.4 in Rybicki & 
Lightman 1979), we find this emission to be in the X-ray band. 



335 keV 



1 - 0.29 



B 



5 X 1013 G 
Mqs\ /lO km 



Mp 







R 



QS 



(7) 



nl/2 



where gravitational redshift is included in the last term. This intensity decays temporally in a 
way closely resembling SGR bursts (see Fig. 6). The decay profiles of our simulated bursts depend 
on the initial plasma /? at the surface of the star and on the dipole inclination 9. Varying these 
parameters in our model allows us to fit the exponential decay shape and duration of observed 
bursts. 



2.3. Periodicity in emission 

The oscillations shown in Fig. 7, can be interpreted as magnetic "breathing" modes. In this 

situation, regions where magnetic reconnection occurs cause a disparity in the magnetic pressure 
between neighbouring regions. This low pressure region will draw plasma in from the surrounding 
regions causing either more reconnection or enhanced density where, in either case, there will be 
oscillations in the emission (i.e. UeB"^). 

Furthermore, since our simulation is performed in the co-rotating frame, our model is too 
simplistic to reproduce a spin-modulated pattern (e.g. the 8.0 s period observed in the March 
5 event) superimposed on a smooth exponential decay in the light curves. However, as can be 
seen in Figs. 2-3, there are hotter regions with more magnetic reconnection events than others, 
and so spin modulation should occur naturally. We note that a few of these hot spots appear 
simultaneously at random locations, implying spin-modulation may consist of even smaller sub- 
pulses, therefore producing many harmonics in the light curve. In other words, if our model is a 
correct representation of SGRs, observations could constrain the number of hot spots. 



3. Subsequent bursts and quiescent phase 

We note that since the CFL phase transition occurs only once in a given star, there is only 
one giant burst in our model. However, if we assume some form of heating, then it is possible that 
the quark star can revert to a non-CFL state wherein superconductivity is lost, and the interior 
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magnetic field lines are no longer constrained to the vortices. Then presumably the dynamics that 
initially caused the magnetic dipole axis to misalign with the rotation axis, we speculate, should 
once again ensue to some extent. Thus we allow for the possibility of subsequent misalignment 
of the magnetic axis that, upon further cooling back below the deconfinement temperature, can 
trigger the same process as before to give subsequent bursts. One might expect these bursts 
strengths to depend on the amount of heating and resulting misalignment. Possible sources of 
subscqTicnt heating can include accretion from a companion, impact from the accretion of a small 
body, or the quark star passing through a higher density region in the ISM. This leads to the 
prediction that radio emission should pick up again slightly before the subsequent bursts, provided 
there is a favourable line-of-sight. The accurate measurement of the temperature in the star in 
these subsequent bursts would now allow observations of the deconfinement temperature, which 
has eluded QCD physicists so far. 

The quiescent phase is due to vortex expulsion from spin-down and subsequent annihilation 
through magnetic reconnection near the surface. The number of vortices decrease slowly with spin- 
down leading to continuous, quiescent energy release which can last for 10^ to 10^ years (see Ouyed 
et al. 2004 for more details). 

4. Model predictions and observational tests 

In the light of the results presented above we will now discuss our model predictions and offer 
our interpretation of the open issues listed in § 1.1: 

• Following reorganization of the outer magnetic field and its alignment with the rotation axis, 
our model naturally predicts the suppression of further radio pulsations. Our model can thus 
in principle explain why SGRs and AXPs stop pulsating following their bursting era. 

• The high-B (magnetar-strength) radio pulsars that show no evidence of enhanced X-ray emis- 
sion, can be accounted for in our model if they are just neutron stars that have not experienced 

the Meissner effect. As long as the object is a quark star in a superconducting phase, mag- 
netic field decay and reorganization will take effect regardless of the field strength. We thus 
predict that SGR-like bursts with moderate-strength magnetic field may be discovered in the 
future. 

• If observations do confirm that the progenitors of SGRs are very massive stars (Gaensler et 
al. 2005) this would strengthen our model since massive progenitors are more likely to lead 
to massive neutron stars, for which it should be easier to reach deconfinement densities in the 
core following accretion or spin-down. 

• The period clustering of observed SGRs and AXPs can be explained in our model as the time 
necessary for the parent neutron star to cool down sufficiently to experience a transition to a 
quark star. 
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Other predictions from our model are: 

• The association of tlic SGR source (tlic quark star) with a parent radio-pulsar. In other words, 
in at least a few cases (if beaming is favorable) a parent radio-pulsar should be detected in 
the same location in the sky as the SGR before the bursting activity. 

• If observations show quiescent emission before and after the burst then this is a subsequent 
burst in our picture, meaning the SGR will have been a quark star for some time. In this 
case the association with a parent neutron star will be less obvious. 

• Assuming synchrotron dominated emission, we should observe a peak at - 280 keV (5/5 x lO^^ 
for a solar mass quark star with a 10 km radius. 

• The total energy in bursts from moderate-strength magnetic field (lO^^G) quark stars would 
be much weaker according to Eq. 6. We predict burst energy from these stars to be of the 
order of lO^'^erg. 

• Finally, we note that during the stars phase transition from hadronic to quark matter, its 
radius changes from 30% - 50% (Ouyed et al. 2002), naturally resulting in a significant 
increase of the magnetic field at the surface. Therefore, a quark star with a strong magnetic 
field is not necessarily the result of a parent star with a strong magnetic field. So, if quark 
stars undergoing the Meissner effect are indeed the cause of SGRs or AXPs, then stronger 
field strengths in these objects should be expected. 

5. Conclusion 

We presented 3-D simulations of the reorganization of the magnetic field surrounding a newly 
born quark star. The reorganization is a consequence of the star entering the CFL phase, confining 
the interior field to vortices, and leaving the exterior field in its tilted dipole configuration. In 
our model the bursting activity in SGRs and AXPs is explained as magnetic reconnection events 
occurring while the exterior field aligns itself with the interior one. One of the obvious consequences 
of the final state (aligned rotator) is the suppression of the radio-emission in SGRs and AXPs 
following their bursting era. We should emphasize that magnetic fields alone cannot be responsible 
for the properties of SGRs and AXPs, but rather a compact star experiencing the Meissner effect 
is required. A further consequence is that while a \ow-B quark star entering the "Meissner phase" 
would burst in X-ray, this will not necessarily be the case for a high-i? neutron star. While quark 
stars have not yet been observed in nature, our model seems to account for many observed features 
in SGRs and AXPs, thus warranting further investigation. 

We thank Ralph Pudritz, Chris Pethick, and Kaya Mori for discussions. B. N. thanks C.I.T.A. 
for hospitality, and Sigma-Xi for its Grant-in-aid of research. The research of R. O. is supported 
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Alberta Ingenuity Fund (AIF). 



A. Magnetic energy of a vertical field in a sphere 

Consider a strictly vertical magnetic field B = Bz{x,y)z in a sphere of radius R. Since the 
radial component of B must be continuous at the surface, the field inside the sphere must be 

B,{r,^,ip)= ' '^^ , Al 

cos fc) 

where Q is the colatitude at the surface, which is related to r, via 

rsiwd = RsinQ = s , (A2) 

with s the cylindrical radius. To calculate the magnetic energy inside the sphere, we integrate B^ 
over z, then transform the s integral in one over B, using Eq. (A2). After a bit of algebra, we find 
the remarkable expression 

2-7r TT ^ 

^mag = R'^ I dip [dQ sine ^^^f'®''^^ . (A3) 





B. Initial magnetic field configuration 

The initial magnetic field represents the phase where the Meissner effect has forced the field 
inside the neutron star to be strictly vertical, but this reorganization has not yet had time to affect 
the external field. 

In order to keep magnetic energy finite, we have to allow for a transition layer between the 
vertical internal field and the inclined dipolar external field. We are thus looking for a magnetic 
field configuration that minimizes total magnetic energy under the following requirements: 

1. For r > i?2, the magnetic field is that of a dipole of dipole moment rh, inclined by with 
respect to the vertical axis 

2. For r < Ri, the field is strictly vertical. 

Minimizing the magnetic energy in the transition layer R\ < r < R2 implies that B is a potential 
field in that region. 

To find the minimum energy field satisfying these requirements, we represent the magnetic 
field through scalar potentials S and T: 

B = -V X (f X VS) - f X Vr . (Bl) 
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Both the vertical field and the inclined dipole field do not have a toroidal part T, thus we set T = 
everywhere. For the poloidal scalar potential S, we make the ansatz 



(B2) 



where ip) are the spherical harmonics of order 1 and degree m. Equation (B2) is compatible 

with both the inner and the outer field; any spherical harmonics of higher order would only increase 
the total energy, so we do not include them. 

Minimizing the sum of the energy (A3) [with R = Ri] and the energy of the potential field 
for i?i < r < i?2, we find the following coefficients for S after a straight-forward, but somewhat 
tedious calculation: 

„m _ ~ am um _ t frm 



where 



and 



Ait 



=F- 



COS1?2 

sin 



V2 



m = 
m = ±1 



(B3) 



(B4) 



b±i 



oo = < 



a±i 



bo = < 



r_ 
i?i 






Rl-Rl 






1 



Rf 



Rl-Rf 



,r < Ri 

,Ri<r <R2 

,r > i?2 
,r<Ri 
,Ri<r<R2 

,r > i?2 

,r < Ri 

,Ri <r < R2 

,r > R2 

,r<Ri 
,Ri<r<R2 
,r > R2 



(B5) 



(B6) 



(B7) 



(B8) 



To obtain the magnetic vector potential (which we need for the Pencil Code), we use the 
formula 

A=-xxVS . (B9) 
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For a thin transition layer of thickness £ = R2—R1 <^ Ri, the contributions to the magnetic 
energy are 

E{ri<r<R2) = ^^^£^2 + 0(1), (Bll) 
E{r > R2) = 2cos2 t?2^2 + 0{e) . (B12) 
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Fig. 2. — "Mock magnetic field lines" in the xz-plane (rotation axis is along z) at t = 1, 25, 100, 
300 going left to right, top to bottom. Note the Alfven front at t = 1, which is accompanied by 

reconnection events. Reconnection events become more apparent with time as the surface magnetic 
field re-adjusts itself. In the last frame we see an aligned dipole which would mark the end of radio 
activity. 




Fig. 3. — "Mock magnetic field fines" in the xy-plane (down the rotation axis) at t = 1, 25, 100, 
300 going left to right, top to bottom. 
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Fig. 4. — Magnetic energy lost over time. The change in magnetic energy (a in Eq. 6) can be 
estimated from these plots to be between 0.4 — 0.6, implying that energies of order 10^^ erg are 
released during the event. 
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Fig. 5. — Total intensity emitted (PqIs) versus time for /3 = 0.3 and 9 = 15°, 30°, 60° (angle between 
rotation and initial magnetic dipole axis). We note that the higher the inclination angle, the longer 
it takes for the outer magnetic field to align itself. This can be seen in the figures where the tail 
takes longer to flatten for higher angles. Furthermore, the high inclination run shows oscillations 
which can be linked to stronger reconnection events. 
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Fig. 6. — Total intensity emitted (PqIs) versus time for 9 = 15° and /3 = 10, 1,0.1. We note that 
the higher /?, the slower the re-adjustment and flattening of the tail. The decrease in variability 
with increasing P can be explained by a decrease in amplitude of magnetic "breathing" modes. 
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Fig. 7. — Time evolution of the Fourier-power spectrum (contours, right axis) for /3 = 1 and 9 = 15° 
along with the normalized intensity (solid line, left axis). The sub-pulses in intensity are shown 
here to be modulated at ~ 3.25 times the spin-period. In an observer's frame, this modulation 
should be further modulated by the spin period. 




Fig. 8.— 



Sketch of magnetic field inside and outside the sphere. 



